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Ferromagnetic resonance (FMR) is used to study magnetic anisotropy of GaMnAs in a series 
of Gai-^Mn^As/Gai-yAlyAs heterostructures modulation-doped by Be. The FMR experiments 
provide a direct measure of cubic and uniaxial magnetic anisotropy fields, and their dependence on 
the doping level. It is found that the increase in doping - in addition to rising the Curie temperature 
of the Gai-iMiijAs layers - also leads to a very significant increase of their uniaxial anisotropy field. 
The FMR measurements further show that the effective <?-factor of Gai_iMn x As is also strongly 
affected by the doping. This in turn provides a direct measure of the contribution from the free hole 
magnetization to the magnetization of Gai-^Mn^As system as a whole. 
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Magnetic anisotropy in ferromagnetic (FM) semicon- 
ductors such as Gai-ajMn^AsiSi^ is expected to play a 
key role in future spin-based devices based on these ma- 
terials. Although it is now well established that the mag- 
netic anisotropy of IIIi-^Mn^V alloys originates from 
the anisotropy of the valence band^*^ the correlation 
between the magnetic anisotropy and the hole concen- 
tration is not yet well established. 

It has recently been found that doping the Gai_„ALAs 
barriers of Gai-^Mn^As/Gai-yAlyAs heterostructures 
by Be acceptors leads to a significant increase of the Curie 
temperature Tc of the Gai-^Mn^As layerii It was also 
shown that ferromagnetic resonance (FMR) can be used 
for directly determining the magnetic anisotropy of thin 
FM films&S In this paper we use FMR to show that mag- 
netic anisotropy in ultra-thin modulation-doped GaM- 
nAs films changes rapidly with doping level; and that the 
effective g-factor of GaMnAs is also strongly affected by 
the doping. This last finding directly reflects the contri- 
bution of the free hole magnetization to FMR dynamics. 

Ferromagnetic Gai-^Mn^As/Gai-yAlyAs het- 
erostructures were grown on semi-insulating (001) 
GaAs substrates by molecular beam epitaxy (MBE), as 
described in detail in Ref. Three heterostructures 
were used in the current study, all three consisting 
of a 5.6 nm Gai-^Mn^As layer (x = 0.062) followed 
by a 13.5 nm Gao.76Alo.24As barrier doped with Be 
starting at the distance of 1 monolayer away from the 
Gai_ a; Mn 2 ;As layer. The Be flux was kept constant 
during the growth, but the thickness of the doped region 
d-Be was varied: ds e = (undoped control Sample 
#1), 5.3 nm (Sample #2) and 13.2 nm (Sample #3). 
The FMR measurements were carried out at 9.38 GHz 
using a Bruker electron paramagnetic resonance (EPR) 
spectrometer. The experimental set-up and the polar 
coordinate system used in the subsequent discussion 
were described in detail in Ref. 0- Each heterostructure 



was cleaved into three 2mm x 2mm square pieces with 
edges along the [110] and [110] directions, and the square 
pieces were mounted in the EPR bridge with the [110], 
[110], or [010] directions pointing vertically. With the 
dc magnetic field H in the horizontal plane, this allowed 
us to map out the FMR for H at any angle Or between 
H||[001] (normal to the layer plane) and three in-plane 
orientations: H|| [110], [110], and [100], following the 
same procedure as in Ref. |9j. 

Because the GaMnAs layers under consideration are 
extremely thin (~ 6 nm), magnetization measurements 
by SQUID were found to be inaccurate. Instead, we made 
use of the fact that the anomalous Hall effect (AHE) is 
dominated by the magnetization M, and can thus serve 
as a measurement of that parameter. To obtain the value 
of M, in our analysis we assumed that AHE is dominated 
by side-jump scattering - i.e., that M oc Kjiall /^sheet> 
where Knaii is the Hall resistance and K s heet is the sheet 
resistance when H is applied perpendicular to the layer. 10 
Typical magneto-transport data R#az; /R-sheet are shown 
in Fig. 1 for several temperatures. Note that at T = 
4.22K the magnetization saturates at fields of about 5 
kOe for undoped material (upper panel), and at higher 
fields (above 6.5 kOe) for modulation-doped samples 
(bottom panel), indicating that the anisotropy field has 
been modified by the doping. 

The inset in Fig. 2 shows FMR spectra at 4. OK 
for a modulation-doped GaMnAs/GaAlAs:Be film (sam- 
ple #3) in four basic configurations: H||[001], H||[110], 
H|| [110], and H||[100]. Strikingly, sharp FMR peaks are 
observed in all configurations (and persist up to Tc), in- 
dicating strong long-range FM coherence of the Mn ++ 
spins. We find this remarkable, since the 5.6nm-thick 
Gao.94Mno.06As film is equivalent to only one monolayer 
of Mn ions randomly distributed over the specimen. As 
shown in Fig. 2, for intermediate orientations of H be- 
tween H|| [100] and H||[001] the FMR peak H R shifts 
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from 1 kOe to 10 kOe. Additionally, an EPR peak is 
observed around 3.3 kOe for all field orientations, orig- 
inating from a small fraction of isolated paramagnetic 
Mn ++ ions with g = 2.00. By their strong dependence 
on crystal geometry, the FMR spectra in Fig. 2 thus es- 
tablish, that magnetic anisotropy plays a major role in 
determining the fields at which the resonances occur' 
The magnetic anisotropy parameters of a GaMnAs film 
can be obtained by analyzing the angular dependence 
of Hr using the following equations and the coordinate 
system defined in Ref. 0. For tpjj — 45° [H and M in the 
(1T0) plane], 

(w/ 7 ) 2 = [H R cos(9 H -e) + (-4nM + H 2± + H 4± /2 
-#4U /4) cos 20 + {H 4± /2 + ff 4 || /4) cos 40] 
x [H R cos{9 h -9) + (-4nM + H 21 _ + H 4 \\/2) cos 2 9 
+ (ff 4 _ L + i? 4 ||/2)cos 4 0-H 4 ||]; (la) 

and for tp H = 0° [H and M in the plane (010)], 

(w/ 7 ) 2 = [Hr cos(9 h -9) + {-AttM + H 2 ± + H 4± /2 
-H 4ll 1 2) cos 26> + (H 4X /2 + II A /2) cos 40] 
x [H R cos(0 ff - 0) + {-AirM + H 21 _ - 2H i \\) cos 2 
+ {H 41 _ + H 4 \\) cos 4 + if 4 ||]- (lb) 

Here H 2 ± and H 4 ± represent, respectively, the uniax- 
ial and the cubic anisotropy fields perpendicular to the 
film; the anisotropy in the film plane is given by the cu- 
bic field H411; uj is the angular microwave frequency; and 
7 = gfiR^ 1 is the gyromagnetic ratio, g being the spec- 
troscopic splitting factor and h the Planck constant. To 
simplify the analysis, we have ignored the small in-plane 
uniaxial anisotropy field H 2 \\ associated with the differ- 
ence between the [110] and [110] axes' 

To determine the parameters appearing in Eq. (1), 
we first analyze the highly-precise values of FMR fields 
Hr for the three high-symmetry directions (H parallel 
to [100], [110], and [001]) £ An independent determina- 
tion of the g-factor and the three anisotropy fields H 2 ± , 
Hm\ and H 4 ± could not be accomplished from the anal- 
ysis of the these values of Hr alone without additional 
constraints, i.e., it was possible to find nearly identical 
fits for several values of g and anisotropy fields within 
experimental error. To reduce the number of fitting pa- 
rameters, we have first imposed the value of g = 2.00 
that corresponds to individual Mn ++ ions. With this 
constraint, the data for Sample #3 at 4 K yield unique 
solutions of H 2 ± = -4319 Oe, H 4 \\ = 739 Oe, and H 41 _ = 
-1933 Oe. Using these values in Eq. (1), we then obtain 
the angular variation of Hr shown by the dashed line in 
Fig. 3. The dashed curves clearly depart from the data, 
indicating that the assumption of g = 2.00 was not valid. 

On the other hand, we note that - due to the large 
compressive - the effect of the cubic H 4 ± term is ex- 
pected to be completely overshadowed by H 2 ± , and may 
be neglected. In our second approach we will therefore 
assume i? 4 j_ = 0, allowing g, iJ 4 || and H 2 ± as fitting 



parameters. With this approach, the data for Sample 
#3 yield g = 1.80, H 4{{ = 720 Oe and H 2 ± = -5887 Oe. 
Using these values (and H 4 ± = 0) in from Eq. (1), an 
excellent fit to the angular variation of Hr is obtained, 
as shown by the solid curve in Fig. 3. We now use 
the above results as starting parameters to carry out a 
weighted nonlinear least squares fit to FMR positions for 
all values of 0#, allowing all four parameters (g, H 4 \\, 
H 2 ±, and H 4 ±) to vary. The results for Sample #3 are: 
g = 1.80 ± 0.02, H 4J _ = 8 ± 110 Oe, H 4 \\ = 735.3 ± 
20 Oe and H 2± = -5764 ± 90 Oe. Not that the relation 
between the three anisotropy, \H 4 ±\ <C I-H411I <C |i?2J_|, 
confirms our assumption that H 4 ± can be neglected as 
a first approximation. Comparing these rigorous results 
with the parameters obtained from Hr observed for the 
three high-symmetry orientations (H parallel to [100], 
[110], and [001]) under the assumption that H 4 ± = 
shows that the two approaches lead to very similar val- 
ues. In our analysis of the data observed as a function of 
temperature we will therefore use the simpler approach. 
Finally, we note that some anisotropy of the g-factor is 
expected in 2D quantum structures'' However, a fit 
obtained by replacing g with g= (g 2 sin 2 + g 2 _cos 2 0) 1 / 2 
cannot be distinguished from the fit with an isotropic g- 
factor. We have therefore accepted an isotropic g-factor 
as an adequate approximation. 

Returning to Fig. 3, it is clear from the inset that this 
contribution of the holes to the g-factor is enhanced as 
hole concentration increases (i.e., the fits departs further 
from the g = 2.00 curves as the doping level increases). 
The fact that all our attempts to fit the data consistently 
lead to an effective g-factor smaller than 2.00 indicates a 
finite contribution of the orbital magnetic moment to the 
magnetization. This can be understood as follows. The 
total magnetization of GaMnAs has two components: a 
contribution from Mn ++ ions, with their pure spin mo- 
ment corresponding to g — 2.00; and the free hole contri- 
bution, which includes both spin and orbital parts. The 
magnetic moments of the hole spins are described by the 
Luttinger parameter k. Since the p-d exchange integral 
(3Nq > and n > are widely accepted for GaMnAs, 5 
the moments of the holes align themselves in the same 
direction as the moments of the Mn ++ ions. In contrast, 
the orbital part of the hole moment (determined by Lan- 
dau diamagnetic currents) create, through the spin-orbit 
interaction, a magnetization opposite to that of Mn ++ 
ions£ Thus the effective g-factor determining the preces- 
sion of the total M is a weighted average of g-factors for 
the Mn-ion (gMn) an d for the hole (gh), described by the 
expression' 

g ef f = (M Mn + M h )/(M Mn /9Mn + M h /g h ). (2) 

Now in the modulation-doped heterostructures, we ob- 
tain a significant increase of the free hole concentration 
compared to "normal" GaMnAs, which automatically en- 
hances their effect X Owing to the antiferromagnetic p-d 
exchange in GaMnAs between Mn ++ and the holes, the 
value of g w 1.80 can only be achieved by assuming that 
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the magnetization of the free holes Mh is opposite to that 
of the Mn ++ ions, and is positive, as indeed predicted 
in Ref. With this constraint, Eq. (2) gives M h w - 
0.15M Mm 9h ~ 5, thus indicating that Landau currents 
make significant contributions to the magnetization of 
the free holes, and must therefore be taken into account 
in III-Mn-As materials. 

The g-factors and related magnetic properties for T = 
10K are listed in Table I for all three specimens under in- 
vestigation. Although the sheet carrier densities p^-Hall 
obtained from Hall measurements are not rigorously valid 
due to the AHE contribution, they nevertheless provide a 
useful indication of the relative level of the doping. 1 - The 
values of the hole concentration obtained from Tc using 
the mean field model,i p v -MF, are also listed in Table I 
for comparison. One should note that the concentration 
of Mn (both substitutional and interstial) is the same in 
all samples, because deposition of the modulation-doped 
GaAlAs:Be layer after the GaMnAs layer does not affect 
the GaMnAs layer which was already in place The 
observed changes listed in Tabe I for all three samples 

- decrease of the <?-factor with doping, enhancement of 
the uniaxial anisotropy field H2±, and reduction of Hm\ 

- can thus only be ascribed to changes in the hole con- 
centration. 

Measurements of FMR up to the Curie temperature 
Tc enable us to determine temperature dependences of 
both the magnetic anisotropy fields and of the <?-factor 
in the modulation doped samples. These quantities, ob- 
tained using the four basic FMR geometries shown in 
the inset in Fig. 2 and assuming H^± — 0, are plotted 
in Fig. 4 for Samples #1 (undoped; open symbols) and 
#3 (modulation doped; solid symbols). As shown in Fig. 
4(a), FMR occurs above the g — 2.00 resonance position 
(horizontal dash-dotted line) when H is perpendicular to 
the film, and below that position for in-plane H orien- 
tations. Shifts from the dash-dotted line gradually de- 
crease - and eventually vanish - as one approaches Tc. 
But clearly the modulation-doped sample has a much 



stronger shift (the strongest shift observed in any GaM- 
nAs samples studied by FMR) than the undoped sample 
when H is normal to the film, indicating a large increase 
of magnetic anisotropy due to the doping. 

Figure 4(b) illustrates several basic features of mag- 
netic anisotropy and its dependence on temperature and 
on the free hole concentration. First, we note that the cu- 
bic anisotropy fields decrease very rapidly with increasing 
T, while i?2_L drops off much more slowly. And second, 
modulation doping unambiguously increases the perpen- 
dicular uniaxial anisotropy field i?2_L , while reducing the 
in-plane cubic field Hm . These observations are consis- 
tent with theoretical calculations predicting changes of 
magnetic anisotropy with hole concentration, 5 although 
at this point the agreement is only qualitative. Finally, 
changes of the ^-factor seen in the inset - although not 
systematic generally - show a clear trend to decrease in 
doped samples below ^60K. While these (/-factor values 
are only approximate, the low-temperature decrease of 
the g-factor in the doped samples may reflect the fact 
that larger numbers of hole spins from the GaAlAs bar- 
rier will couple with Mn ++ spins in GaMnAs as T de- 
creases, thus increasing the effect of Landau diamag- 
netism on the overall magnetization. 

In summary, the results reported in this paper clearly 
point to the crucial role which holes play in determin- 
ing the magnetic anisotropy and the magnetization of 
Gai-^Mn^As. It is especially important to note that 
magnetic anisotropy of IIIi-^Mn-r V materials can be ma- 
nipulated by controlling the hole concentration, thus pro- 
viding a mechanism which may be employed in future 
device applications. Since the dependence of the mag- 
netic anisotropy parameters on the free hole concentra- 
tion is not yet well understood, further rigorous theoret- 
ical studies addressing this issue are clearly needed. 

This work was supported by the DARPA SpinS Pro- 
gram, by the Director, Office of Science, Office of Ba- 
sic Energy Sciences, and by NSF-NIRT Grant DMR02- 
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TABLE I: Key parameters for the Cai^Mn^ As/Gai-y Al M As 
heterostructures modulation-doped by Be. 



Sample # 


1 


2 


3 


Structure 


GaMn As / GaAlAs 


GaMnAs/GaAlAs:Be 


dse (nm) 





5.3 


13.2 


Tc (K) 


72 


85 


95 


p v -MF (cm- 3 ) 


1.24xl0 20 


1.48 xlO 20 


1.64xl0 20 


p s -Hall (cm -2 ) 


1.32xl0 14 


1.74xl0 14 


2.94xl0 14 


M WK (emu/cm 3 ) 


30.0 


25.9 


29.4 


HIT (Oe) 


-3734 


-5560 


-5924 


HIT (Oe) 


827 


471 


557 


gl0 K 


1.94 


1.87 


1.82 




Applied Magnetic Field [kOe] 



FIG. 1: Magnetization curves deduced from AHE at various 
temperatures for Samples #1 and #3. Magnetic field H is 
applied along the hard axis of magnetization, H|| [001] . 




dc Magnetic Field H (kOe) 



FIG. 2: FMR spectra (T = 4K) observed for Sample #3 at 
various orientations 8h (from 0° to 90° in 10° increments) for 
H between [100] and [001] directions in the (010) plane. The 
dotted line is a guide for eyes, indicating the shifting FMR 
position. The insert shows FMR spectra for the perpendicular 
(H||[001]) and three parallel (H||[110], H||[lT0], and H||[100]) 
configurations observed at 4.0 K. The weak peaks observed 
at the g = 2.00 resonance position (indicated by the vertical 
dashed line) are ascribed to EPR of isolated paramagnetic 
Mn ++ ions. 
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FIG. 3: Angular dependence of FMR positions for Sample #3 
for H in the (110) plane (right-hand panel), and for H in the 
(010) plane (left-hand panel). Inset: Angular dependence of 
FMR positions at 4.0 K for H in the (lTO) plane for all three 
samples used in this study. Dashed curves show theoretical 
fits obtained for g = 2.00, H$± ^ 0. The solid curves are fits 
obtained for Hi±_ = 0, g = 1.80 (main figure); and Ha± — 0, 
g =1.80, 1.87, and 1.92 (inset, top to bottom). 
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FIG. 4: Temperature dependence of the FMR results for un- 
doped (Sample #1, open symbols) and modulation-doped 
(Sample #3, solid symbols) GaMnAs/GaAlAs heterostruc- 
tures. Figure 4(a) shows FMR positions observed for the four 
basic orientations of H (same as in the inset in Fig. 2). Figure 
4(b) shows uniaxial anisotropy fields i?2± for the two samples 
(bottom panel); cubic anisotropy fields H411 (top panel); and 
the effective g-factors (inset). 



